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Abstract 

A blueshift tunable meta materia I (MM) exhibiting a double-negative refractive index based on a topological 
insulator (bismuth selenide, Bi 2 Se 3 ) has been demonstrated in the near-infrared (NIR) spectral region. The potential 
of Bi 2 Se 3 as a dielectric interlayer of the multilayer MM is explored. The optical response of elliptical nanohole arrays 
penetrating through Au/Bi 2 Se 3 /Au films is numerically investigated using the finite difference time domain (FDTD) 
method. The blueshift tuning range of the MM is as high as 370 nm (from 2,140 to 1,770 nm) after switching the 
Bi 2 Se 3 between its trigonal and orthorhombic states. 
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Background 

Metamaterials (MMs) are artificially engineered compos- 
ites that attract considerable interests due to their excep- 
tional electromagnetic properties, which are not typically 
found in nature, such as negative refractive index and 
cloaking [1-4]. These MMs with various subwavelength 
resonant elements have offered magnetic and/or electric 
resonant responses to incident electromagnetic radia- 
tion, scalable from the microwave frequencies up to the 
terahertz and optical ones [5-7]. Particularly, nanohole 
resonators embedded in metal-dielectric-metal (MDM) 
multilayers are frequently used as building blocks of 
negative- refractive- index MMs [8-11], owing to the cou- 
pling between surface plasmons counterpropagating on 
the two closely spaced interfaces which results in a 
closed loop of the electric currents. This gives rise to 
magnetic dipolar resonances between the two coupled 
metal layers, while the continuous metallic strip parts 
provide the electric resonance moments [12,13]. All these 
features make the nanohole array perforating through 
MDM films become a strong candidate for developing 
three-dimensional negative-index MMs [14,15]. 

One of the obstacles in this progress is the resonance 
responses of MMs to the impinge light which are usually 
fixed once the dimension of the structure is determined, 
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thus making the MMs possess a limited bandwidth. How- 
ever, for many applications (switching, modulation, filte- 
ring, etc.), it would be highly desirable to tune the MM 
resonances over a wide bandwidth. To this end, tunable 
photonic MMs, the spectral range of which can be con- 
trolled by changing the dielectric environment of the res- 
onator with liquid crystals (LCs) [16-18]; phase transition 
materials [19,20]; and optical pumping [21,22] have been 
discussed recently. However, the challenge is to develop 
tunable MDM-MMs in the near-infrared (NIR) regime. It 
is due to the fact that frequency tunability of the MDM- 
MM mainly requires for the interlayer dielectric material 
to possess a tunable effective dielectric constant in the NIR 
region, hence limiting the choice of the active materials. 
Here, we take a different approach to actively tune the res- 
onant frequency of the MDM-MMs in the NIR regions by 
using bismuth selenide (Bi 2 Se 3 ) as the dielectric layer. 

Recently, a rising Dirac material - topological insulators 
(TIs) - had been intensively researched in condensed mat- 
ter physics [23,24]. In analogy to the optoelectronic 
applications of graphene, a thin layer of TIs has been 
theoretically predicted to be a promising material for 
broadband and high-performance optoelectronic devices 
such as photodetectors, terahertz lasers, waveguides, and 
transparent electrodes [25]. Among these TIs, Bi 2 Se 3 is a 
particularly interesting compound due to its relatively 
large bulk band gap and a simple surface state consisting 
of a single Dirac cone-like structure [26,27]. Study of 
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the dielectric function reveals that the optical dielectric 
constant of Bi 2 Se 3 can be very different for the trigonal 
and orthorhombic phases in the NIR regime [28]. Bi 2 Se 3 
exhibits a number of means through which their dielectric 
properties can be altered [28-33]. Herein, structural phase 
transition between trigonal and orthorhombic states, 
which is achieved by a high pressure and temperature, is 
proposed and studied as a means to change the intrinsic 
effective dielectric properties of the MDM-MMs [28] . 

Here, we numerically demonstrate a blueshift tunable 
nanometer-scale MM consisting of an elliptical nanohole 
array (ENA) embedded in the MDM multilayers where 
the dielectric core layer is a Bi 2 Se 3 composite. Under a 
high pressure of 2 to 4.3 Pa at 500°C, Bi 2 Se 3 occurring 
in trigonal phase undergoes a transition to orthorhombic 
phase and features a large change in the values of the effec- 
tive dielectric constant [28]. Accordingly, a massive blue- 
shift of the resonant response (from 2,140 to 1,770 nm) of 
a Bi 2 Se 3 -based MDM-ENA is achieved in the NIR region. 
Our proposed blueshift tunable negative-index MM pro- 
vides greater flexibility in the practical application and has 
a potential of enabling efficient switches and modulators 
in the NIR region. 

Methods 

The proposed MDM-ENA suspended in a vacuum is 
shown in Figure 1, with the coordinate axes and the pola- 
rization configuration of the normally incident light. The 
structure consists of trilayers of Au/Bi 2 Se 3 /Au. The thick- 
ness of each Au layer is 30 nm, and the thickness of the 
Bi 2 Se 3 layer is 60 nm. The metamaterial parameters are op- 
timized for the maximum sensitivity of the resonance to a 
change in the refractive index of the Bi 2 Se 3 core dielectric 
layer in the NIR spectral range. The element resonator is 



shown in Figure lb, where the pitch of the elliptical holes 
is L = 400 nm, the diameters of the elliptical holes are 
di = 240 nm and d 2 = 120 nm, and ft is a cross-sectional 
plane of the structure. The z-axis is normal to the struc- 
ture surface, and the x-y plane is parallel to the structure 
surface. This simulated structure is periodically extended 
along the x and y axes. The tunable optical properties of 
the structure are calculated using 3D EM Explorer Studio 
[34], a commercial finite difference time domain (FDTD) 
code. In the simulation, a simple Drude-type model for Au 
permittivity was used, which is a good approximation to 
experimental values in the NIR region. The dielectric prop- 
erties of Au as given by Johnson and Christy are used [35]. 
A plane wave source is simulated at normal incidence to 
the structure. The computational domain (400 nm x 
400 nm x 1,000 nm) has a perfectly matched layer (PML), 
absorbing boundaries in the z direction and periodic boun- 
daries in the x-y plane [36]. A uniform FDTD mesh size is 
adopted. The mesh size is the same along all Cartesian 
axes: Ax = Ay - Az = 2 nm, which is sufficient to minimize 
the numerical errors arising from the FDTD method. 

The topological insulator material Bi 2 Se 3 was selected 
due to its significantly different optical properties be- 
tween the trigonal and orthorhombic phases. The real 
(£i) and imaginary (s 2 ) parts of the dielectric function 
for the different structural phases of Bi 2 Se 3 were ob- 
tained from the published data in [28]; the NIR spectral 
region is shown in Figure 2. A large change in the di- 
electric function across the NIR is obtained after swit- 
ching Bi 2 Se 3 from trigonal to its orthorhombic phase. 

After the complex coefficients of transmission t — T a 
e l< Pa and reflection r = R a e l<Pm are obtained by the 3D EM 
Explorer Studio, in which T a is the amplitude and (p a is 
the phase of the transmission coefficient, and R a is the 




=120nm 



(a) (b) 

Figure 1 Schematic of the proposed structure, (a) Schematic of the MDM structure consisting of a 60-nm-thick Bi 2 Se 3 dielectric layer between 
two 30-nm-thick Au films perforated with a square array of elliptical holes suspended in air. The lattice constant is L = 400 nm, and hole diameters 
are d-\ = 240 nm and d 2 = 1 20 nm. (b) Illustration of the square lattice of ENA. 
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Figure 2 Dielectric constant of the Bi 2 Se 3 . (a) Real part of dielectric function e-i (oj) for trigonal and orthorhombic phases of Bi 2 Se 3 . (b) Imaginary 
part of dielectric function e 2 {d) for trigonal and orthorhombic phases of Bi 2 Se 3 . 



amplitude and (p ra is the phase of the reflection coeffi- 
cient, the effective optical parameters can be extracted 
using the Fresnel formula [37]. 

For an equivalent isotropic homogenous slab of thick- 
ness h surrounded by semi-infinite media with refractive 
index n x and n 3 under normal incidence, we have 
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The so-called material parameters £ e ff and ^eff of a sin- 
gle layer of a double-fishnet negative-index metamaterial 
are extracted using the well-known Nicholson-Ross- 
Weir (NRW) method [38-40]. Therefore, once n eff and rj 
are evaluated, the effective permittivity and permeability 
are calculated using 



(3) 



where n e & is the effective refractive index, rj is the im- 
pedance, h is the thickness of the structure, k = o/c, c 
is the speed of light, yh is an arbitrary integer, and vi\ — 
n 3 = l since the structure is suspended in a vacuum. The 
signs of n e ff and rj and the value of m are resolved by the 
passivity of the metamaterial that requires the signs of 
the real part of impedance rj and imaginary part of ef- 
fective index # e ff to be positive, i.e., Real(^) > 0, Imag 
(« ef f ) > 0 which is consistent with the study described in 
[39,40] . We then apply this extraction approach to deter- 
mine the change in the optical response of the structure 
when the phase of Bi 2 Se 3 is switched between its trigonal 
and orthorhombic states. 

Results and discussion 

The ENA has a lower transmittance for s-polarized light 
due to the electric field's orientation with respect to the 
metallic stripe width [12]; hence, the polarization of the 



incident wave was set to be j?-polarized. As shown in 
Figure la, s polarization means that the incident electric 
field vector is parallel to the long axis of the ENA, and 
the incident electric field vector perpendicular to the 
long axis of the ENA is then denoted by p polarization. 
We first investigate the transmittance T= \ t\ 2 and re- 
flectance R=\r\ 2 of the structure for p polarization in 
Figure 3. Structures with a different dielectric constant 
of Bi 2 Se 3 (shown in Figure 2) were modeled to investi- 
gate the effect of the phase change of Bi 2 Se 3 on the 
position and amplitude of the spectrums. It can be 
seen that the resonance wavelength blueshifts from 2,140 
to 1,770 nm when the structural phase of Bi 2 Se 3 swit- 
ches from trigonal to orthorhombic. The structure is 
impedance-matched, hence possessing a low reflectance 
corresponding to the dips in reflectance of Figure 3b 
for different forms of Bi 2 Se 3 . 

In Figure 4, the transmission (t) and reflection(r) pha- 
ses are demonstrated. The transmission phase exhibits a 
dip around the resonance, indicating that the light is 
advanced in phase at the resonance, characteristic of 
a left-handed material [41], Importantly, changing the 
structural phase of the Bi 2 Se 3 offers transmission and 
reflection phase tunability which implies tunable ef- 
fective constitutive parameters in the structure. 

Taking into account the subwavelength thickness of 
the structure, the extracted n ei f can be retrieved from 
the transmission and reflection coefficients shown in 
Figure 5. For the MM with the trigonal Bi 2 Se 3 dielectric 
layer, the negative-index band extends from 1,880 to 
2,420 nm with a minimum value of the real part of the 
refractive index Real(# e ff) = -7- Regarding losses, the 



RealQeff) 



is ta- 



figure of merit (FOM) defined as FOM — Imag ( w ff ) 
ken to show the overall performance of the MM, where 
Imag(w e ff) is the imaginary part of the refractive index. 
As shown in Figure 5c, the FOM for the trigonal phase 
is 2.7 at the operating wavelength of 2,080 nm. The 
negative- index band of the orthorhombic Bi 2 Se 3 -based 
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Figure 3 Transmittance and reflectance. 3D FDTD simulation of (a) spectrum of transmittance and (b) spectrum of reflectance, for the 
different phases of the Bi 2 Se 3 dielectric layer, where the light source is p polarization at normal incidence angle. 



MM extends from 1,600 to 2,214 nm having a minimum 
value of Real(^eff) = -3.2. The FOM is 1.2 at the reson- 
ant wavelength of 1,756 nm. Furthermore, the band- 
width of Real(^eff) becomes wider for the orthorhombic 
Bi 2 Se 3 film in Figure 5a due to increased damping of the 
plasmon resonance but at the cost of an accompanied 
lower value of FOM in Figure 5c. 

The refractive index is expressed in terms of the real 
and imaginary parts of the permeability ^ e ff and per- 
mittivity £ e ff- However, the sign of the real part of the 
permeability ^ eff : Real(^ eff ) determines the relative mag- 
nitudes of the imaginary and real parts of the refractive 
index [41]. To achieve a negative index with a small loss, 
a negative Real(^ e ff) is required. Therefore, we have sim- 
ulated ^ e ff and £ e ff for the structure as shown in Figure 6. 
For the trigonal and orthorhombic phases of Bi 2 Se 3 , Real 
(^ ef f) has a Fano-type line shape and Im(^ eff ) has a 
Lorentzian line shape in the region of the negative index. 
Moreover, a double-negative MM can be achieved when 
Real(^ e ff) and Real(£ eff ) simultaneously reach negative 
values over a wide frequency range and thus a reduced 
loss. The maximum negative Real(^ e ff) decreases with 
the phase transition from the trigonal to orthorhombic, 



hence resulting in the smaller value of the maximum 
negative Real(w ef f) at the orthorhombic phase. 

This magnetic negative response can be explained loo- 
king at the current and field distribution at the reso- 
nance wavelengths. Figure 7 shows the current and total 

magnetic field intensity H = y\H x \ + \H y \ + \H Z \ for 

the magnetic resonant wavelengths of 2,140 and 1,770 
nm at the /? plane shown in Figure 1. In the field maps 
of Figure 7, the arrows show the currents, whereas the 
color shows the intensity of the magnetic field. It clearly 
shows that the antiparallel currents are excited at op- 
posite internal metallic interfaces, closed by an electric 
displacement current / D . Therefore, these virtual current 
loops between two Au layers on the ft plane give rise to 
magnetic resonant responses of negative Re^eff) that 
interact strongly with the incident magnetic field at 
which the total magnetic field intensity H is strongly 
localized in the Bi 2 Se 3 dielectric layer between the top and 
bottom Au layers [42] . 

Specifically, H for the orthorhombic phase shown in 
Figure 7b is weaker than the trigonal phase shown in 
Figure 7a. It depicts that the MM based on orthorhombic 
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Figure 4 Transmission and reflection phase. 3D FDTD simulation of (a) phase of transmission and (b) phase of reflection, for the different 
phases of the Bi 2 Se 3 dielectric layer, where the light source is p polarization at normal incidence angle. 
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Figure 5 Effective index and figure of merit. 3D FDTD simulation of (a) real part of n effl (b) imaginary part of n effl and (c) figure of merit for 
the different phases of the Bi 2 Se 3 dielectric layer, where the light source is p polarization at normal incidence angle. 
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Figure 6 Permittivity and permeability. 3D FDTD simulation of (a) the real part of /L/ e ff, (b) the imaginary part of /L/ e ff, (c) the real part of £ e ff# 
and (d) the imaginary part of £ eff for the different phases of the Bi 2 Se 3 dielectric layer, where the light source is p polarization at normal 
incidence angle. 
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Figure 7 Magnetic field intensity and displacement current. A map of the normalized total magnetic field intensity distribution (color bar) 

and displacement current (J D ) (arrows) along the plane /3: (a) at a 2,140-nm resonance wavelength for trigonal Bi 2 Se 3 and (b) at a 1,770-nm 

resonance wavelength for orthorhombic Bi 2 Se 3; where the light source is p polarization at normal incidence angle. 
^ J 



phase has a smaller magnetic dipolar moment than the 
trigonal phase and thus smaller FOM. 

To further understand the negative -index resonance in 
the metamaterials, it is useful to study the dispersion of 
the surface plasmon polariton (SPP) modes within the 
multilayer structure. Both the internal and external SPP 
modes in the multilayer metamaterials are similar to 



those of the same structure without resonant elements, 
i.e., MDM films [42], where the internal SPP mode reso- 
nates in the inner surfaces of the metal layers and the 
external SPP mode resonates in the outer surfaces of the 
metal layers. Therefore, the SPP dispersion relation of 
the multilayer metamaterial can be approximately ap- 
proached by that of the MDM structure. In Figure 8, we 
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Figure 8 Dispersion relation of the structure. Representation of the dispersion relation of the Au-Bi 2 Se 3 -Au trilayer (left) and the transmittance 
of the multilayer metamaterials (right) for both (a) trigonal Bi 2 Se 3 and (b) orthorhombic Bi 2 Se 3 . 
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have calculated the SPP mode dispersion relation of the 
Au-Bi 2 Se 3 -Au sheets with the top Au film thickness 
ti - 30 nm, middle Bi 2 Se 3 film thickness t 2 = 60 nm, 
and bottom Au film thickness t 3 - 30 nm. The trans- 
mittance spectrum of the multilayer metamaterials is 
also depicted together with the dispersion relation of 
the Au-Bi 2 Se 3 -Au films. 

Recalling the coupling condition from light to SPP 
modes [42], it can be seen that the (1,1) internal reson- 
ance of the Au-Bi 2 Se 3 -Au trilayer is excited at 2,350 nm 
associated with the trigonal Bi 2 Se 3 in Figure 8a. This in- 
ternal SPP resonance blueshifts to 2,010 nm when the 
trigonal state changes to the orthorhombic state as shown 
in Figure 8b. We also observe that the two internal (1,1) 
modes which appear at 2,350 and 2,010 nm in the simple 
MDM structure do not perfectly match the two absorb- 
ance peaks at the resonance wavelengths of 2,140 and 
1,770 nm in the multilayer metamaterials for both the tri- 
gonal and orthorhombic phases, respectively. This differ- 
ence is because the dispersion relation of the SPP modes 
used as matching condition does not include the resonant 
squares, which cause a resonance shift [42]. 

Conclusions 

In conclusion, this work numerically demonstrates the 
tunable optical properties of an ENA perforated through 
Au/Bi 2 Se 3 /Au trilayers. We present that the MDM-ENA 
can be improved to exhibit a substantial frequency tun- 
ability of the intrinsic resonance in the NIR spectral 
region by selecting Bi 2 Se 3 as the active dielectric ma- 
terial. Particularly, the resonant transmission, reflec- 
tion, and the effective constitutive parameters of the 
Bi 2 Se 3 -coupled multilayer MM can be massively blue- 
shifted by transiting the phase of the Bi 2 Se 3 film from 
the trigonal to orthorhombic. This may offer an in- 
novative and practical paradigm for the development 
of tunable photonic devices. We expect that our re- 
sults will facilitate further experimental studies of the 
tunable MMs and make this technique suitable for 
tuning of plasmon resonance in the optical regime. 
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